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Abstract 

We present explicit expressions for the Mellin transforms of Laguerre and 
Hermite functions in terms of a variety of special functions. We show that 
many of the properties of the resulting functions, including functional equa- 
tions and reciprocity laws, are direct consequences of transformation formulae 
of hypergeometric functions. Interest in these results is reinforced by the fact 
that polynomial or other factors of the Mellin transforms have zeros only on 
the critical line Re s = 1/2. We additionally present a simple-zero Proposition 
for the Mellin transform of the wavefunction of the D-dimensional hydrogenic 
atom. These results are of interest to several areas including quantum mechan- 
ics and analytic number theory. 
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Introduction 

The work of Bump et al. El [TT] on a property of the zeros of the Mellin trans- 
forms of Hermite and associated Laguerre functions has created significant interest. 
Since the zeros of these functions lie on the critical line Re s = 1/2, suggestions of 
connections with the Riemann hypothesis have been made. The Mellin transforms 
of orthogonal polynomials have additional properties and some of these have been 
pointed out [5, 6J. 

In this paper a complementary point of view is developed through the theory 
of special functions. We obtain explicit results with the aid of, and in terms of, a 
variety of special functions. In particular, when the Mellin transforms are written 
in terms of the Gauss hypergeometric function, well known transformation formulae 
yield functional equations, reciprocity laws, and other properties. 

Although many of our final results recover those found in Refs. [5] and [6], our 
approach may be useful in finding clues to other algebraic connections, since, for 
instance, the harmonic oscillator Hamiltonian may be given a group theoretic in- 
terpretation via the Weil representation of SL^R) [TT]. Our work is related to the 
development of extended theta function representations of the Riemann zeta function 
£. We have very recently shown how to correspondingly generalize the important 
Riemann-Siegel integral formula [9]. Bump and Ng [5] and Keating p3] have shown 
how to generalize Riemann's second proof of the functional equation of £(s) by using 
Mellin and Fourier transforms of Hermite polynomials. These types of results extend 
to L-functions and automorphic forms [6l [TJJ . As we encounter the classical Whit- 
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taker function, there are analogs in a more general automorphic context. Much of 
this appears to have an interpretation in terms of group representations. 

We recall that the eigensolution of the fundamental quantum mechanical prob- 
lems of the harmonic oscillator and hydrogenic atoms contain Hermite or associated 
Laguerre polynomials, depending upon the coordinate system used and the spatial 
dimension (e.g., [2Tj [8]). Two quantum numbers appear for indexing the energy 
levels and the angular momentum. These wavefunctions have a wide variety of ap- 
plicability, including to image processing and the combinatorics of zero- dimensional 
quantum field theory [8]. Very recently additional analytic properties of these "quan- 
tum shapelets" have been expounded [8]. The one- dimensional Coulomb problem has 
recently reappeared as a model in quantum computing with electrons on liquid he- 
lium films [231122] . In addition, given the self-reciprocal Fourier transform property of 
Hermite polynomials, there are several applications in Fourier optics [7J [13] . Hermite 
and Laguerre polynomials are also important in random matrix theory. 

The simple harmonic oscillator and Coulomb problems may be transformed to one 
another. For example, the 4-dimensional (8-dimensional) harmonic oscillator may be 
transformed to a 3-dimensional (5-dimensional) Coulomb problem [T6l IT8] . Many 
more mappings are realizable, especially for the corresponding radial problems |16j . 
Therefore, the eigensolutions are closely related for these two problems with central 
potentials. 

We specifically call attention to Ref. [3] wherein Mellin transforms were calculated 
for hypergeometric orthogonal polynomials and relations within the Askey scheme [15] 
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were discussed. However, that reference did not consider the location of the complex 
zeros of the transforms or connections to the integral representations of zeta func- 
tions. Moreover, given a family of orthogonal polynomials p n (x, {«j}) where are 
real parameters, orthogonal with respect to a weight function w(x, {oti}), we are most 
interested in the Mellin transform of the associated functions yjw(x, {ai})p n (x, {ctj}). 
This manner of including the weight function is very useful in identifying the orthog- 
onality and other properties of the Mellin transforms, and was often not the case in 
Ref. [3]. 

We first consider the Mellin transformation of Laguerre functions, outlining mul- 
tiple proofs of a result in terms of the hypergeometric function 2 F\. We then briefly 
present similar results for the Mellin transforms of Hermite functions. Corollaries 
of our Propositions include the functional equations and reciprocity laws at negative 
integers of polynomial or other factors of the transforms. Other results include gener- 
ating functions and recursion relations of the Mellin transforms. We also strengthen 
a zero result of Ref. [6] . We give a simple- zero result for the Mellin transform of the 
wavefunction for .D-dimensional hydrogenic atoms and then supply some concluding 
remarks. 

Mellin transform of associated Laguerre functions 

We put for the Laguerre functions C%(x) = x a ' 2 e~ x / 2 L%(x) for a > —1, where L" 
is the associated Laguerre polynomial (2J [12J, [19] . We also put for the Mellin transform 

roc 

M°(s) = / C a Jx)x s ~ l dx = 2 s+a ' 2 T{s + a/2)P«{s). (1) 
Jo 
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This transform exists for Re (s + a/2) > 0. 
Proposition 1 We have 

p:{s) = <±±°k 2Fl (-n, s + a/2; a + 1; 2), (2) 

where (a) n = T(a + n)/T(a) is the Pochhammer symbol and T is the Gamma function. 
Corollary 1 The polynomials -P°(s) satisfy the functional equation 

TO = (-irai-4 (3) 

Corollary 2 The polynomials -P"(s) satisfy the reciprocity law 

^^P-i-m - a/2) = ( -^^P«(-n - a/2). (4) 
ml n\ 

Corollary 3 The polynomials P"( s ) have derivative 

Corollary 4 The polynomials -P"(s) may be written as 

= t TT^y U-^ +1 s(k,j)(s + a/2)', (6) 

where s(k,j) are Stirling numbers of the first kind. 

Equation (3) has the form £(s) = a£(l — s) where \a\ = 1. Such a functional 
equation is standard in the theory of completed zeta functions. Corollary 1 follows 
from Proposition 1 upon the use of the transformation formula 



2 F 1 (a, (3; j; z) = (1 - 2 F X (a, 7 - # 7; . (7) 



Corollary 2, that has a combinatorial interpretation from counting lattice points 
[6], follows from Proposition 1 from the obvious symmetry 2Fx(—n, —m; a + 1; 2) = 
2-Fi(— m, —n; a + 1; 2). Corollary 3 follows from Proposition 1 by using the series def- 
inition of 2-Fi and the derivative of the Pochhammer symbol in terms of the digamma 
function if) = T'/T. We have 

^-P:(s) = { -±±^ ± -t^f (s + a) k ^{s + a/2 + k) - ^(s + a/2)]. (8) 
as n\ fc=Q [i -+- a)k K\ 

Applying the functional equation of the digamma function [T2] then gives Eq. (5). 
The same method may be used to obtain higher order derivatives of P^{s). The series 
representation (6) for -P"(s) follows from Proposition 1 by the series definition of 2F1 
and the expression for a Pochhammer symbol in terms of s(k,j) pQ. 

In fact the polynomials of Eq. (2) are closely connected with the symmetric 
Meixner-Pollaczek polynomials P^(x,ti/2) [31 [15]. The latter polynomials are a 
special case of P^ x \x, <p) where A > and < <p < ir. For the polynomials of Eq. (2) 
we identify P£(s) = (-i) n P^ l+a ^ 2 \i/2 - is,n/2). 

For Proposition 1 we describe six separate proofs. These are in addition to the 
equivalent series result of Ref. [6] presumably based upon the power series of the 
associated Laguerre polynomials [T2l [T9] and termwise integration. We believe that 
the alternative proofs may contain illuminating intermediate results that may have 
other applications. Additional proofs are possible using, for example, other integral 
representations of L". As this work was being finished, we found that Eqs. (1) and 
(2) are effectively given in Ref. [20] (p. 245). Our alternative proofs illustrate a 
variety of analytic techniques and exhibit connections with other special functions of 
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mathematical physics, in line with the theme of this paper. 

First proof. We start with the definition of M* and make a change of variable: 

roc 

M%(s) = 2 a ' 2+s / e-v y a l 2+s - l L a n {2 y )dy. (9) 
J o 

We then apply the property \±0\ 

I&rx) = £ ( P + m )r n (l-rr~-L^x) (10) 
t^o \ m ~ n J 

at t = 2, giving 

n ( en A- n\ r 00 

M;(s) = 2 a / 2+s (-l)" ]T(-1)^' _ / e- y y a/2+s - 1 LJ{y)dy. (11) 

The integral on the right side of this equation exists for Re (s + a/2) > and is given 
in terms of [12] (p. 844) 2-Pi( — j, a/2 + s; a + 1; 1). In turn, this function value may be 
reduced to the ratio of Gamma functions [12J (p. 1042) by Chu-Vandermonde sum- 
mation. Making these substitutions, transforming the Gamma factors, and writing 
them in terms of Pochhammer symbols, gives Eq. (1) with P£(s) as stated in Eq. 
(2). 

Second proof. Here we start with an integral representation of the associated 
Laguerre polynomials in terms of Bessel functions of the first kind J a and develop an 
intermediate integral representation of M n (s) in terms of a Whittaker function M M) „. 
We have [2] (p. 286) 

P x /2 poo 

£«( x ) = __^ t n+a/2 J a (2v / ^)e-^, (12) 
so that the interchange of integrations gives 

1 roo roo 

M?(s) = — / t n+a l 2 e~ l dt / x s - 1 e x/2 J a (2V^i)dx. (13) 
n\ Jo Jo 
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Since pj (p. 720) 

T(s + a/2),f iy /2 - 



/ x s - l e x ' 2 J a {2Vrt)dx = A_t L2 e '(--) M s _ 1/2 , a/2 (-2t), (14) 
Jo vtr(ot + l) V 2/ 



with a change of variable we have 



km = r f;, + ^ r ^ +a/2 - i/2 M,_ 1/2 , a/2 (y)^ (is) 



n!r(a + 1) 7o 

= r ,w + ^ (-l) n+Q/2 ~ s+1 2 a/2+ T(n+a+l) 2 F 1 (n+a+l, l-s+a/2; a+1; 2). (16) 
n!l (a + 1) 

In the last step we used Ref. [T2J (p. 859). Examining the asymptotic form of M^ u (z) 
for large \z\ shows that formula 7.621.1 there is valid for Re s > —1/2, not just for Re 
s > 1/2. The reduction of Eq. (16) to Eq. (1) with (2) is accomplished by applying 
the transformation formula 



2 F 1 (a, (3- 7; z) = (1 - zf-^ 2 F l { 1 - a, 7 - # 75 4 (17) 
Third proof. From the power series of the Laguerre polynomials we have 

MZ(s) = (1 + a) n £ ,,, 4^T- , / ^/ 2+s - 1+fe e-^ 2 dx. (18) 
^ fc!(n - fc)!(l + a) k Jo 

We evaluate the integral in terms of the Gamma function. We then manipulate the 

series, using l/(n— k)\ = (—n) k (—l) k /n\, and Eq. (2) follows. 

Fourth proof. We substitute into the definition of M"(s) the expansion [19] (p. 

89) 

= r( , + a /2)r ( , - a /2) g r(j - - 1 + ^ L L_ _ _ , (19) 

make a change of variable, and twice apply property (10). We then use the orthogo- 
nality of the associated Laguerre polynomials with respect to the measure x a e~ x dx. 
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Fifth proof. We make use of the connection with the confluent hypergeometric 
function X F X (e.g., [19J, p. 273), 

L «( x ) = (l±^k 1 j? 1 (_ n;a;a .). (20) 

We then use a tabulated integral [12] (p. 860), giving Eq. (2). Thus, this method is 
among the most expedient in obtaining the desired Mellin transform. 

Sixth proof By using the exponential generating function of L° [2] (p. 288), it is 
readily verified that [6] 

oo 

Y,M%(s)t n = 2 s+a/2 T(s + a/2)(l-t) s - a/2 -\l + t)- s - a/2 , \t\ < 1. (21) 

n=0 

By then using the generalized binomial expansion 

00 (a) 

(l-f)- = £^f» = ^iM;M), (22) 

n=0 

we have 

oo 

P%{s)t n = iF (a/2 - s + 1; .; t) iF (a/2 + s; .; -t) (23a) 



n=0 

oo 



= E E — ~ + S) '"' (-l)"-'f. (23.) 

n=0£=0 *-V" 

By comparing powers of t on both sides of Eq. (23b) we have 

^ £!(n-£)! 

(-l) n r(«/2 + n + s) . . . , ,, 

2 Fi(-n, l + a/2-s;l - a/2 - n - s; -1). (246) 



n!r(s + a/2) 

We then apply the linear transformation [20J (p. 48) 

2 iq(a, 6, c; z) = ^, x~/*~l> ~ a 2^1 (a, a - c + 1; a + 6 - c + 1; 1 - - 
T{c — a)T{c — b) V -2 
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r(c)T(a + b-c) 



+ v ^(1 - zf- a - b 2 Fi ( c - a, 1 - a; c - a - b + 1; 1 - 



(25) 



r(a)r(6) 

The appearance of the prefactor 1/T(— n) before the second resulting 2 F t an d the 
existence of simple poles of the Gamma function at nonpositive integers annuls the 
second term. We then simplify the first term using T(—a) /Y(—a — n) = (1 + a) n and 
Eq. (2) follows. 

We remark that Eqs. (24) may be obtained as a highly degenerate case of a 
formula 2 F 1 (a,b; c; gz) 2 Fx(a, (3; 7; hz) = Y^= c k zk i where the coefficients may be 
expressed in terms of a generalized hypergeometric function 4 F 3 . In our case, b = c, 
(3 — 7, and g = —h = 1. 

We next present example recursion relations satisfied by M"(s) and -P"(s). We 
have 

Proposition 2. 

1, 



: [M«( S ) + M« +1 ( S )] = (a/2 + s - 1)K +1 ( S - 1) - M-( S - 1)], 



i 1 

1 + - - + s 



n V2 



1 



M«(s) = —M%(s + !)+! + - M^s 



2n 



a 



n 



(s + a/2)[P:(s) + P: +1 (s)] = (a/2 + a - l)[P n Q +1 (s - 1) - P n a (s - 1)], 



(26a) 

(266) 
(27a) 



and 









1 ~i 




(HI 




n 





P«( S ) = -(a + a/2 + l)P:(s + !)+! + - P^s). (276) 



n 



a 



77 



Proof. Equations (27) follow from Eqs. (26) upon using the definition of P"( s ) i n 
Eq. 1. Equation (26a) follows from the use of the relation [12] (p. 1037) 

d 



dx 



[LZ(x) - L« n+1 (x)\ 



(21 
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in Eq. (1) and integration by parts. Similarly, Eq. (26b) follows from the relation 
[32] (p. 1037) 

^Ki^M + C 1 * < 29 > 

and integration by parts. 

Bump et al. [6] noted the orthogonality of the polynomials P°(l/2 + it) with 
respect to the measure 2 a+1 \T(l/2 + a/2 + it)\ 2 dt. By applying the Plancherel formula 
they demonstrated that all the zeros of -P"(s) have real part 1/2 ([6], Theorem 4). 
By applying classical results [25] (Theorems 3.3.1-3.3.3), we may strengthen their 
conclusion to 

Proposition 3. (a) The zeros of P%(s) are simple and lie on the line Re s = 1/2. 
(b) The zeros of P°(l/2 + it) and P" +1 (l/2 + it) separate each other. 

Proposition 3 also follows from Theorems 5.4.1 and 5.4.2 of Ref. [2]. 

Due to the relations Z^ 1/2 (x) = (-l) n H 2n (^) /2 2n n\ axidL]l 2 (x) = {-l) n H 2n+l {y/x) 
/y/x2 2n+1 n\ [19J (p. 81), where H m are the Hermite polynomials, the theory of the 
Hermite polynomials may be deduced from that of the associated Laguerre polyno- 
mial in the cases that a = ±1/2. Nonetheless, we believe that some separate 
discussion for Hermite polynomials is in order. We develop the hypergeometric func- 
tion representation of their Mellin transforms, relate these to other special functions, 
and consider reciprocity relations for a factor of the Mellin transforms. 
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Mellin transform of Hermite functions 

Following the normalization of Ref . [5] , we put 

f n (x) = (8n)- n / 2 H n (V2^x)e-™\ (30) 

For the Mellin transform, we set 

M n (s) = 2 / f n (x)x s - l dx = n- s / 2 T(s/2)p n/2 (s), Re s > 0. (31) 
Jo 

We first consider the case of even degree Hermite polynomials and have 
Proposition 4. 

Pn (s) = (8*)-»(-l)»M: 2 F 1 (-n, S /2; 1/2; 2). (32) 

Corollary 5 The polynomials p n (s) satisfy the functional equation 

Pn (s) = (-l) n Pn (l-s). (33) 

Corollary 6 The polynomials Pn (s) satisfy the reciprocity law 

(8n)- m (-ir^p n (-2m) = (87r)^(-l) n ^p m (-2n). (34) 
ml n\ 

Corollary 7 The polynomials p n {s) may be expressed as 

Pn(s) = (87r)- ^^ (2n)!C|/ 2 -(V2), (35) 

where is the Gegenbauer (generalized Legendre) polynomial [121 [13 ED] • 

Corollary 5 again follows from the transformation (7). Corollary 6, that may not 
have been observed before, follows from the symmetry of 2 Fi when its two numerator 
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parameters are interchanged. It would be of interest to know if Eq. (34) has a 
combinatorial interpretation. Corollary 7 follows from the relation [12] (p. 1030) 



where B is the Beta function. In this way, known recursion formulas for the Gegen- 
bauer polynomials may be used to obtain the corresponding ones for the polynomials 

Pn(s). 

A great many methods may be used to obtain Eq. (32). We simply indicate two 
of them here. 

Proof 1. If we use the relation pj] (p. 1033), H 2n (x) = (-l) n (2n)! iFi(-n, 1/2; x 2 )/n\, 
with a change of variable we may use a tabulated integral p2] (p. 860). 
Proof 2. If we perform the Mellin transform directly, we have 



= 2(87r)- n (27r)- s/2 ^^2 M/2 - (s - 1)/2 4 s/2 r(s/2)r(n + 1/2) 2 F 1 (-n, s/2; 1/2; 2), 



where we used [12] (p. 838). Equation (32) is obtained after simplification and noting 
that r(n + 1/2) = ^{2n)\2~ 2n /n\. 

The function p n /2{s) is no longer a polynomial when n is an odd integer. We now 
have 

Proposition 5. 



Pn+1/2 = 2^r n - 1,2 {-mn + l/2)M r(g ^ + y 2) 2 Fx(-n, s/2 + 1/2; 3/2; 2). 

n\ T(s/2) 

(38) 



2 F 1 (-n, s/2; 1/2; 2) = (-l)^B( s /2 - n, n + 1)0*^(^2), 



(36) 




(37) 
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In this case, the reciprocity law relates p n+ i/ 2 (— 2m — 1) and p m+ i/ 2 (— 2n — 1). The 
Mellin transform M2 n+ i(s) may be performed directly, with the aid of [12] (p. 838). 
While p„ + i/ 2 (s) is not a polynomial in s, the truncating 2 F± in Eq. (38) is. 

Bump and Ng [5] (p. 197) gave a recursion relation for Pk{s). Many more recursion 
formulae exist for Pk{s) and M n (s) and we mention two such. We have 
Proposition 6 

1 (n — 1) 

M n (s) = — M n _ x (s + 1) - ^— ^M n _ 2 ( S ), (39) 

V27T 47T 



and 

MJs 



(s - 1) 

M n+1 {s + 1) - ^—^M n+1 (s - 1) 

Z7T 



(40) 



4(n+ 1) 

Equation (39) follows from the use of H n (u) = 2uH n -i(u) — 2{n — l)if n _ 2 (?i) in Eq. 
(1) while Eq. (40) follows from the use of H n (u) = H' n+1 {u)/2(n + 1) [T2] (p. 1033). 

Finally, we obtain the generating function of M n (s) and p n /2(s) and describe how 
the functional equation (33) arises from that of the parabolic cylinder function D v . 
By using the generating function of the Hermite polynomials [19J (p. 60) we have 

MJs) 

n=0 

2(87r)- n/2 e {n - 1)t2 (2Tt)- s/2 T(s)D^ s (-2t), Re s > 0, (41) 



Z^n Til JO 



where we used [12] (p. 337) to evaluate the integral. Thereupon, from Eq. (31), we 
have 

£ P^ditr = 2(87r)-"/ 2 e^- 1 ) i2 ^Pr ( — ) D_ 5 (-2t), (42) 



n=0 711 ^ 
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where the duplication formula of the Gamma function has been used. The important 
functional equation (33) is recovered from this equation by using that of -D_ s [20] 
(p. 325). Another way to see this is to note the connection with the confluent 
hypergeometric function [20] (p. 324) 



D„{z) = 2^V 22/4 



1(1 / 1/2; z 2 /2) + i^[(l - i/)/2; 3/2; z 2 /2) 



(43) 



and that Kummer's first transformation ii r i(a; p; z) = \F\{p — a; p; —2) applies. 

Mellin transform of the solution for hydrogenic atoms in D-dimensions 

We now return to some of the quantum mechanical considerations of the Introduc- 
tion. We let ip(x) be the wavefunction for the hydrogenic atom (Coulomb problem) 
in D- dimensions, where x = (x\, . . . , xd) and r = \x\. We then have 
Proposition 7. The Mellin transform 

tfj(x)r s - D/2 ~ 1 dx (44) 



IRP 

has zeros only on the critical line Re s = 1/2 and these zeros are simple. 

Proof. For hydrogenic atoms of nuclear charge Ze, the scaled Hamiltonian is given 
by H = —V 2 + V(r) with potential energy V(r) = —Ze 2 /r, where e is the electronic 
charge. The eigensolutions satisfy Hip ni = E n ip n i, where n > 1 is the principal 
quantum number, £ — 0, 1, . . . ,n — 1 the angular momentum quantum number, and 
E n oc —l/rj 2 (n) the energy levels, with 77(72) = n + (D — 3)/2. (The energies are 
degenerate, meaning that they are independent of I here.) For the central potential, 
the wavefunction is separable, so that ip(x) = YR n i(r), where the function Y is 
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independent of r, and explicitly we have [21] 



77(71) 

where the normalizing constant K n p is independent of r. We have 



(45) 



P P POO 

/ ip(x)r q dx = / YdQ / RneMrir^dr, (46) 
Jr d Js - 1 Jo 

where S 1 ™" 1 is the sphere in R n and dQ the measure on it. We put Eq. (45) into 
the radial integration of this equation and change variable to u = 2r/r](n). We apply 
Proposition 3 with q = s — 1 — D/2 and Proposition 7 follows. 

Remarks, (i) We have just considered the standard Coulomb problem in D di- 
mensions. However, it should be noted that only in three dimensions is the potential 
V(r) oc 1/r the same as the Greens function of the Poisson equation — V 2 = 47rp, 
where <fi is the electrostatic potential and p the charge density [21]. (ii) The rela- 
tivistic pi-mesic atom is described by the Klein-Gordon equation. This problem may 
be transformed to the hydrogenic atom in any dimension, and the wave equation 
solutions are functionally identical [21]. In particular, the same radial wavefunc- 
tion (45) appears. Therefore, Proposition 7 could instead be stated in terms of the 
wavefunction of a pi-mesic atom, (iii) The 3-dimensional hydrogenic atom is also 
separable in parabolic coordinates (£,77,0) [Ml Hj- In this case the wavefunction 
u nin2m (Z,ri,<P) oc e-«/ 2 e |m|/2 47U(0e-^ H/2 i:S m (r/)e^, m = 0, ±1, ±2,..., 
contains Laguerre functions in both of the confocal paraboloid coordinates. 
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Summary and brief discussion 

We have explicitly evaluated the Mellin transforms of Hermite and associated 
Laguerre functions in terms of several other special functions. In particular, our 
approach permits the arsenal of known results for hyper geometric functions to be 
applied. For instance, the transformation formulae and symmetries of the function 
2F1 directly lead to the functional equation and reciprocity law of the polynomial 
factor of the Mellin transform. The Hermite and associated Laguerre functions and 
their Mellin transforms are of much interest in analytic number theory due to their 
connections with generalizations of Riemann's second proof of the functional equation 
of his zeta function j5j[6l[T7]. This interest is reinforced since the polynomial factors 
of the Mellin transforms have simple zeros that occur only on the critical line Re 
s = 1/2. 

We have described links of the Mellin transforms to fundamental problems of quan- 
tum mechanics, including the isotropic harmonic oscillator and Coulomb problems. 
Our Proposition 7 evidences such connections for the solution of the Schrodinger equa- 
tion with central potential. Classically, only for the harmonic oscillator and Kepler 
problem are the orbits always closed. Transformation between these two problems 
in quantum mechanics was noted by Schrodinger himself and there is continuing in- 
terest in this subject. In addition to our work, known algebraic relations (e.g., [IB] ) 
may yield additional insight into the group representation aspects behind the Mellin 
transforms of interest. We also mention that character sum analogs over finite fields 
are known for the Hermite polynomials [TT] . 
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The orthogonality of factors of Mellin transforms of solutions of the Schrodinger 
equation may be viewed as follows. Since the Schrodinger equation with suitable 
boundary conditions is self-adjoint, its eigenfunctions corresponding to distinct eigen- 
values (eigenenergies) are orthogonal. Then the isometry of the Fourier (Mellin) 
transform takes the orthogonality in real space to that in momentum space. 

Our Proposition 1 suggests several questions. The most begging question may be 
the following. If we put f(s) = 2-^1 s + a/2; a> + 1; 2), Eq. (7) shows its functional 
equation to be f(s) = (—1)^/(1 — s). Therefore, we ask what conditions must be 
placed upon the parameter f3 in order for / to have zeros only on the critical line? 
The development of an orthogonality relation for such functions with respect to an 
appropriate measure would probably address this. At least two other hypergeometric 
extensions may be possible. We first consider all cases of the extended hypergeometric 
function p F q that reduce to polynomials (e.g., [20], p. 64) and correspond to a self- 
adjoint boundary value problem. We then ask which of these have a functional 
equation relating argument z to 1 — z and possess zeros only along Re z — 1/2. 
Another direction is to consider multivariate hypergeometric functions. This includes 
the two-variable Kampe de Feriet function and its generalizations. 

Another avenue of generalization would be through the connection of the extended 
hypergeometric function to the Meijer G- function ([12J, p. 1071). In particular, 
interesting candidate functions are offered by the G- function at argument 1 or —1, as 
then the functional relations of G show there are functional relations at parameters 
a r , b s and 1 — b s and 1 — a r . 
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Within the Askey scheme |15j . the Wilson polynomials W n (x 2 ;a,b,c,d) with Re 
(a,b,c,d) > are orthogonal on [0, oo) with respect to the measure (27r) _1 |r(a + 
ix)T(b + ix)T(c + ix)T(d + ix) /T(2ix)\ 2 dx. Therefore it seems worthwhile to consider 
in the future the Mellin transform of the function 



W n (x 2 ; a, b, c, d) 



r(a + ix)T(b + ix)T(c + ix)T(d + ix) 



V(2ix) 



W n (x 2 ; a, b, c, d). 



(47) 
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